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Abstract: We show spurious effects in perturbative calculations due to different orderings 
of inhomogenous terms while computing corrections to Green functions for two different 
metrics. These effects are not carried over to physically measurable quantities like the 
renormalized value of the vacuum expectation value of the stress-energy tensor. 
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Introduction 

Finding logarithmic behaviour of correlation functions instead of pure power law has always 
been interesting in quantum field theory. QCD differs from a free theory at the asymptotic 
region only by its logarithmic corrections. The presence of these corrections differentiate 
a physically interesting theory from a trivial model. Another effect is the possibility that 
logarithmic corrections in perturbation theory may sum up to an anamolous power law, as 
seen in the Thirring model [1]. Keeping these examples in mind, it is always get exciting 
to encounter logarithmic behaviour while calculating Green functions. That is why one 
may be curious whether encountering such terms always points to important physical 
phenomena. 

Here we want to study the interaction of a gravitational field described by a certain 
metric with a scalar field. We do not know how to quantize the gravitational field, though. 
This fact necessitates the use of semi-classical methods . These methods treat the gravita- 
tional field classically and couples it to the scalar field only by writing the d'Alembertian in 
the background of that metric [2] . This operation reduces the full field theoretical problem 
to an external field calculation. We study the n-point functions of the problem and try 
to deduce information about the full theory from these quantities. Out of these n-point 
functions, the two point function, which is the Green function of the d'Alembertian, is the 
most important one, since by differentiating it we can get the vacuum expectation value, 
VEV, of the stress-energy tensor, T^ v . 

To be able to calculate renormalized value of the VEV of the stress-energy tensor, 
T^u >renj we have to find an algorithm to regulate the divergences of the two-point 
function, Gf, at the coincidence limit. The initial step in this programme is to study 
the singularity structure of Gp for the particular model studied. We know that the sin- 
gularity at the coincidence limit of the Green function for a free scalar field in flat space 
is quadratic. At this point one is confronted with a very important theorem [3] concern- 
ing the singularity behaviour of Green functions in different backgrounds. This theorem 
states that if a metric is flat at any region in space-time, the singularities of the Green 
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function in the background of this metric should exhibit Hadamard behaviour [4], which 
is at worst quadratic. Although we may by-pass the results of this theorem by studying 
metrics that are not C°°, still it is very improbable for Green functions of scalar particles 
in the background metrics that are flat at any region in space-time to have singularities 
that are worse than quadratic. 

The theorem quoted above makes us suspicious of any extra logarithmic terms we 
encounter while calculating the Green functions in the background of different metrics, 
especially if the calculation is performed perturbatively. We want to see if such terms 
are genuine or the result of the different regularization procedures used, arising when we 
make our expressions unnecessarily singular. Whether such effects are genuine or not can 
be checked by comparing perturbative results with the exact ones, when available, and 
studying different ways of grouping terms with more or less singular behaviour. 

Here we are going to continue our investigation of spurious effects in Green function 
calculations using different metrics [5]. We will see that while we take alternative routes 
in solving inhomogenous partial differential equations results in solutions with different 
singularity structure for the Green's function, this difference cancels out in the calculation 
of physical quantities like the VEV of the stress-energy tensor, after it is renormalized. 

Our model is an impulsive gravitational plane wave solution of Hogan [6]. 

4-OT-J ffj J 

ds 2 = -2 cos 2 aU\dz + Q(U) H(z)dz\ 2 + 2dUdV, 1 

(X 

with the only one non-vanishing Ricci tensor component, 

R44 = -2a 2 Q(U), 2 

and only one non- vanishing Weyl tensor component , 

#4 = H{z)8{U). 3 

Here © is the Heavyside unit step function, 5 is the Dirac delta function, a is a constant 
and H is the z derivative of a smooth arbitrary function of z, H(z) = ^£L. U and V are 
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null coordinates, z = -^(x + iy) . This metric is flat for U < and conformally flat for 
U > 0. This metric is similar to the metric found by Nutku earlier [7] . 

We have two reasons for studying this model . First, we have the exact result to 
compare with the perturbative one, in one of the two cases studied. Second, we know that 
the < >ren should be null Z 8 ' 9 , a result we have to obtain at the end. 

We see that this metric has a dimensional parameter a. We have seen that the presence 
of a dimensional parameter in the metric may result in non-standard ( non-Hadamard) 
behaviour in the Green function , when this quantity is calculated perturbatively / 10 . 
Here we find that this change in the singularity structure occurs if we group our terms in 
a certain manner in the perturbation expansion even in the first order calculation. It is 
absent for the exact calculation . It is amusing, however, to see that these terms are absent 
if we group the inhomogenous terms in a different manner. This fact shows that the change 
in the singularity in the former calculation is due to the wrong choice of the parameters, 
which results in severer singularities, and, when regulated in the Schwinger formalism, 
ends up with terms which are worse than the other case by logarithms. We find that this 
is another method for generating logarithmic singularities in the coincidence limit for the 
Green functions at will . A similar phenomenon was shown to exist in [10], by taking the 
homogenous solutions into account. These singularities do not survive when the VEV of 
the renormalized stress-energy tensor, a measurable physical quantity, is calculated. 

We first solve the Green function, Gf for a special H(z) used in the metric given 
by eq. [1] exactly and point to the Hamadard behaviour of this expression. Then we get 
the perturbative solution in two ways, and show two different results obtained for Gf, 
only one of them in the Hadamard form. In the third section we use a more complicated 
form of H(z) , where we could not obtain the exact result and show the same conflicting 
behaviour when Gf is calculated perturbatively. At the end we note that when the VEV of 
the stress-energy tensor is computed the different expressions for G f give the same result 
for < T^ u > ren . All through our work we use conformal coupling. 
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2.1 Exact Calculation for the First Metric 



We choose G = az for the arbitrary function in the metric given above, which gives the 
following expression for the d'Alembertian operator written in the background of this 
metric 

□ = 2-|— - 2atan(2aC7) J- + — J— - (sm(2aU) - 1) °~ 



dUdV v ' dV cos 2 (2aU) y v ' ' dx 2 

4 d 2 
(sin(2a*7) + 1) 



If we write 



we find 



cos 2 (2aU) dy 2 
nf(x,y,U,V) = 0, 



p iRV p ikxx p ik 2 y k*-k% k 2 +k 2 
f - P iRa l os (Lu) P \mr tan(2aC/) fi 

1 cosV2( 2a [/) ' 
Here k±, k2, R are the Fourier modes that will be integrated over in the Green's function 
calculation. These integrations are performed easily and we get 

a(0(U-U>)-0(U>-U)) 



8V2i7v 2 sin 2a(U - U') ({V - V) - 



(x-x') 2 (y-y) 2 



where 



-L— + tan(2aC7) + ] - tan(2a£/')) , 8 

cos(2a(7) cos(2a(7') / 

A i=( 77, — 77V tan(2a£7) -j— — - tan(2aC/') J . 9 

1 V cos ( 2at/ ) cos(2at/') v 7 

We see easily that in the coincidence limit this function has a quadratic divergence, same 

as in the flat metric. 



2.2 Perturbative Calculation for the first metric 

We expand the operator □ in powers of a. If we write 

□ = Lq + aL x + 10 
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we find 

d 2 d 2 d 2 

L ° ~ 2 dUdV ^dx^^ dy 2 "' 11 

L - 2 ( ^ &2 \ 12 

1 \dx 2 dy 2 J 

We also expand both the eigenfunction and the eigenvalue of the equation 

□0 = 13 
in powers of a as <fi = <fi° + a<p l + ... and A = A + A 1 + .... We find 

Ap = 1 c i(AU+RV+kix+k 2 y) 14 

(2n) 2 (2\R\) 1 2 

\° = K-kl-kl 15 

A 1 = (0°, Li4>°) = 0. 16 

Here k\ , hi , K , R are the different modes we have to integrate over to find the Green 
function Gp- satisfies the equation 

(L o -Ao)0 1 +L 1 0° = O. 17 

The form of this equation suggests the ansatz 

<f> 1 = <f>°(g(U,z) + h(U,z)) 18 

where z = -^(x + iy). This ansatz yields the following equation for the unknown functions 
g and h: 

iR {% + w)- > + *>i - - + (kl + kl)U = °- 19 

We seperate this equation into two parts, one for a function of z and J7, and the other a 
function of z and U . One choice for this decomposition is taking 

iR w - 7^ {kl + lk2) ^i + \ {kl + kl)u = °' 20 
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and 



iR^r - ^=(h - ik 2 )^ + \{k 2 + k 2 )U = 0. 21 



At this point a clearification is in order. Since the inhomogeneous term is a constant, as 
far as z and z are concerned, there is ambiguity how it is shared among the two equations. 
Here we designate by gi,h\, the particular choice for the decomposition given above. 

Once this seperation is made, the integration is immediate. We find 

gi = - -M± Uz _ iR (M + g) Jg2> 22 



iV2(k 1 +ik 2 ) ' 2(ki + ik 2 ) 2 ' 



h\ = Wi±MLuz - tR n f 1 + k R 9 z 2 . 23 
y/2(k 1 -ik 2 ) 2(k 1 -ik 2 ) 2 

To calculate the Green function Gp we have to sum over all the modes, 

A 

which reduces in our first order calculation to 

c (i) _ Mo [9i + 9l + h 1 + h* 1 ] 25 

When written explicitly, we get 

/oo pco POO PCO 

dR dK dh / dk 2 x 
-oo J — CO J — CO J — CO 

Ji[R(V-V')+^(U-U')+k 1 (x-x')+k 2 (y-y')]) 

x (2,r mm -k- } -k-, ™ + « 



where 

and x is the generic name for the four variables U, V, x, y. We use the Schwinger represen- 
tation to write 



ik\ k 2 ) f . R 2 ;2w..2 



1 Z" 00 

— ro ro = W (ictexpf— i(K — k 2 — k )a — ai(5), 

K-k\-k\ J 
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, , , , , . exp (-i(Atf + ki)P - 7/?) 29 

where 5 and 7 are infinitesimal positive and real quantities. We perform the integrals in 
the usual manner. The end result is written using an infrared parameter m and the zeroth 
order Hankel function as 

G? = ( ^ + w*. W (WS) + J) , 



4^(2^)2 \BA B 2 
where ... contains terms with the same singularity behaviour as the first two, i.e. terms 

with quadratic divergence and quadratic divergence times a logarithmic divergence in the 
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coincidence limit. The logarithmic divergence is given by the Hq ' term which goes to a 
logarithm as its argument goes to zero when the infrared cut-off is removed, 

H {2) oc log A + log 2m, 31 

and we discard the log 2m term. In the above expression 

A= (x — x'){Ux - U'x') -(y- y')(Uy - U'y'), 32 

B=(x-x') 2 + (y-y') 2 , 33 

A = 2(U - U')(V - V) -{x- x') 2 -(y- y'f . 34 

In this expression we see that the expected Hadamard behaviour, i.e. the worse divergence 
being only quadratic, is modified by a logarithm. 

It is amusing to note that this behaviour, which is not reflected to the exact solution, 
is an artefact of the choice we used in separating our eq. [19] into the holomorphic and 
antiholomorphic parts. Another choice is separating eq. [19] as 

tR W - T2 { 1 + 2) ^7 + 1 u = °' 35 

lR du ~ " lk2) ~m + — I — u = a 36 

The integrations of these equations give immediately 



h 2 = — 2 =(k 1 -ik 2 )Uz-iR^-. 38 
2^2 4 

Note that with the latter choice for decomposition of the equation, we have reduced the 
powers of (k\+ik 2 ), (ki — ik 2 ) in the denominator. Since we integrate these expressions from 
minus infinity to plus infinity, terms in the denominator vanish in this range and we use 
the Schwinger prescription given above to regulate them. The calculation of Gf with less 
severe divergences is much simpler now. The Green function integration is straightforward. 
For U > U', it reads 

Gf = — f^- x 

4^(2^)2 



x 



C 1 [(x - x'){Ux - U'x>) - (y - y')(Uy - U'y')] + C 2 [(u - u')(x 2 - y 2 - x 2 + y 2 )] 



A 2 

39 

where C\ , C 2 are two constants. This singularity has the same singularity behaviour as 
the exact solution at the coincidence limit. 



3. Second metric 

Here we show the same thing with a different function G(z) used to specify the metric in 

eq. [1] explicitly, to illustrate that the phenomena we find is not special to only one choice 

2 

of the trial function H . We take the next simplest form, G = 9 ^-. Then the metric reads, 
in the region where it is not flat, 



ds 2 = 2dUdV - 2 
We find 



dzdz (cos 2 (aU) + a 2 zzsm 2 (aU)) + %in(2aU) (z(dz) 2 + z(dz) 2 ) 
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_ _2acos(a*7)sin(a*7)(l + zz) d d 2 2acos 3 (aU) sm(aU) ( d ; d 

D " b W + dudv + 53 Vai + ai 

2acos(aC/)sin(a(7) / d 2 _d 2 \ 2a 2 cos 2 (aU) sin 2 (aU) ( 8 _d 

Z-p—r + Z-— — r Z— h Z-— 



B 2 V dz 2 dz J B 2 \ dz dz 

^cos 2 (aU) + a 2 zz sin 2 (aU) d 2 ^ 
B 3 dzdz 
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Here B = cos 2 (all) — a 2 zz sin 2 (all). We could not solve the Green function of this operator 
exactly; so, we do not know the exact singularity structure of it. We compare, however, 
the singularity structure of the two expressions we obtain by grouping the inhomogeneous 
terms differently. Just as in the first example, one way results in a function with Hadamard 
behaviour, the other gives rise to a term which is modified by logarithmic corrections. 



We expand the operator □ in powers of a. At the first nontrivial order we get 



□ 



(- 



d 2 



\ dUdV dzdz 



d 2 \ 



+ 2a 2 U 



d d d 
dV dz dz 



+ z 



dr 



+ z 



dz 2 
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We separate □ into two parts □ L + a 2 Li and expand both the eigenvalue A and the 
eigenfunction of the equation 

0(f) = Xcf) 43 

in the same manner. The zeroth-order and first-order eigenvalues and zeroth-order eigen- 
function are as given in eq.s (14)- (16). We get the equivalent of eq. [17] , with new 
L\ for 0i. The structure of this equation again suggests the ansatz of eq. [18], (j>\ = 
4>o[g(U, z) + h(U,z)] which gives rise to the equation 

d (ik 1 -k 2 \dg (ik 1 + k 2 \dh 

lR du {9 +h) ~ {—^r ) dz - {—jr ) & 



+u 



, , ikl Z k2 ]\z( ^t^ V - V2i kl - iR 



= 0. 
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V2 J \ V2 

Here we get two different results depending on how we separate this equation into two 
equations for the two unknown functions. One choice is to write 



. j-,dgi ( %k\ — k 2 \ dgi i Uz r . t , n2 ^-R^ , iUk\ 



iR 



dU 



iR- 



V2 

%k\ + k 2 \ dhi Uz, . 



+ 



V2 



dU \ ^2 J dz 
These equations are integrated to get 



, . . , iRU , iUk\ 

— + —(ih + k 2 ) - - 



+ 



V2 



= 0, 



= 0. 
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iRz 3 



9i = 



+ 



2(k 2 + k 2 ) 2 + 2(k 2 + k 2 ) 2 ~T U ^ k2 > 
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U(iR-iV2k 1 )(ik 1 + k 2 ) 

z 47 



k^ ~~\~ k<2 

and the corresponding expression for h\ where the z is replaced by ~z and ik\ + k 2 goes into 
ik\ — k 2 and vice versa. If we use this solution to obtain the Green function, we obtain 
functions with quadratic singularity at the coincidence point as well as functions whose 
singularities are modified by a logarithmic term. It is also amusing that the logarithmic 
behaviour comes out only when we have an ambiguity in seperating the equation. The 
first , second , fourth and the fifth terms have only the quadratic Hadamard singularity, 
whereas the third and the sixth terms , in this form, give rise to logarithms. The details 
of this calculation can be found in reference [11]. We want to state only that when we 
calculate the Green function for the term that reads 

9 , (ih + fcg) ^ 2 

plus the ~z part , we get 

„ ( %f? ( r~ x L 2 (* 2 + *' 2 -y 2 + y 2 + + 2*y)ff< 2) (2mA) 49 

{{x - x'Y + (y - y'YY 

plus terms with the same singularity structure in addition to terms with only quadratic 
singularity structure. A in this expression was defined in eq. [34]. We can modify the 
third and the sixth term in eq. [47] to reduce the power of the terms in the denominator, 
though, and write the differential equations as 

dg 2 (k 1 + ik 2 )dg 2 Uz 2 RU U 

R du - ~v^~^ ~ ^ {kl + lk2) - — - 71 ( 1 + 2) = °' 50 

dh 2 (fci - ik 2 ) dh 2 Uz 2 RU U 

R mJ - —y^^ " ^ {kl " lk2) - — - 71 ( 1 " 1 2) = °- 5 

In this expression we keep most of the terms same as those given in eq.[45] and change 
only one term , the only term which gave the logarithmic correction. 

The solution of the above equations reads 
iRz 3 . 2 / 'U(k!+ik 2 ) R 2 V2iR \ TT / iR 



6 ' \ V2 2(h+ik 2 ) 2 2(h+ik 2 )J V V2(h + ik 2 ) 

52 
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and the similar expression for h 2 (U,z). We check the behaviour in the coincidence for 
parts of Gp that are replaced in the new expression. The terms that are changed are the 
fourth and the fifth terms in eq. [52] . We can show that these new terms give rise to only 
to quadratic divergence in Gf .The calculation of the fifth term gives rise to 

({x-x^+(y-y'r)A ( ( * " ^ + " ^ " ^ + % " V ' ] { * V + * V) ) 53 

with no logarithm. The calculation of the fourth term gives only the A term in the 
denominator and exhibits Hadamard behaviour . 



4. Conclusion 

We have already noted [10] that while solving eq. [17] after making the ansatz 4>i = 
(j>o(g(z, U) + h(z, U)), we are solving an equation of the type 

rD d (ki + ik 2 ) d 

[R W vT~ & ,9 = / 64 

where / denotes the inhomogeneous part of the equation. This equation also has a solution 
for the case when / = 0, the homogeneous case. Indeed, any arbitrary function with the 
argument + ik 2 ) + y/2z is a solution. If the function / has a different dimension 

as compared to 4>q, then we have to multiply the solution to the homogenous equation 
by a dimensional constant. In our problem the only such quantities are the modes, R 
and ki + ik 2 that exist in eq. [14]. This seems to be completely innocent, as far as g is 
concerned. The fact that we have to sum over all the modes, however, changes the result 

(2) 

obtained for Gp. One shows that we can generate Hq , the function that modifies the 
Hadamard behaviour in this way. Which one of the two factors, R or k\ +ik 2 , is used does 
not change the character of the new singularity structure. It was shown in [10] that only 
Hq survives the two derivatives ifwe want to extract the vacuum expectation value of the 
stress-energy tensor, < T Mi , > from Gf using the established methods / 2 . Itis also shown 
that this new singularity structure which already exists inthe homogeneous solution, and 
which is here unsurfaced by taking different combinations for the inhomogenous term, does 
not change the result that plane waves do not polarize the vacuum / 8 ' 9 . 
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Here we show another way to generate such spurious effects. We refer to our earlier 
work / 10 to establish that the renormalized value of the VEV of the stress-energy tensor is 
independent of all these spurious effects. The essence of the argument is that there is no 
finite part of the resulting expression when the coincidence limit is taken, since we have 
only an isolated pole divergence, and no remaining finite part. 

We see that there is more than one way to generate spurious singularities for per- 
turbatively computed Green functions. Here we study examples which are obtained by 
grouping terms in different manners in inhomogenous differential equations. It is a relief 
that these spurious effects do not change the value calculated for < T Mi , > ren . 

Acknowledgement: We are grateful to Prof. Dr. Yavuz Nutku for giving us his metrics 
and general support. This work is partially supported by TUBITAK, the Scientific and 
Technical Research Council of Turkey. M.H.'s work is also partially supported by TUBA, 
the Academy of Sciences of Turkey. 



13 



REFERENCES 



1. K.G.Wilson, Phys. Rev D2, (1970) 1473 and 1478 ; 

2 N.D.Birrell and P.C.W. Davies, Quantum Fields in Curved Space Space, Cambridge 
University Press, Cambridge (1982) , S.A. Fulling, Aspects of Quantum Field Theory 
in Curved Space-Time Cambridge University Press, Cambridge, 1989, R.M. Wald, 
Quantum Field Theory in Curved Spacetime and Black Hole Thermodynamics, The 
University of Chicago Press, Chicago, 1994; 

3. S.A.Fulhng, M.Sweeney and R.Wald, Comm. Math. Phys. 63 (1978) 257; 

4. J.Hadamard,Lectiires on Cauchy's Problem in Linear Partial Differential Equations, 
Yale University Press, New Haven (1923) ; 

5. Earlier work can be found in M.Hortagsu and R. Kaya , Gen. Rel. Grav. 30 (1998) 
1645; 

6. P.A.Hogan, 1997 " Plane gravitational Waves and Bertotti- Robinson Space-Times" , 
University College Dublin , preprint; 

7. Y.Nutku in Directions in General Relativity, Papers in honor of Charles Misner ed. 
by B.L.Hu, M.P.Ryan Jr. and C.V. Vishveshara, Cambridge University Press , Cam- 
bridge (1993); 

8. S.Deser, J. Phys. A: Math. Gen. 81 (1975) 972; 

9. G.W.Gibbons, Comm. Math. Phys. 45 (1975) 191; 

10. M.Hortagsu and K.Ulker, , Class. Quantum Grav., 15 (1998) 1415; 

11. B.C.Liitfuoglu, I.T.U. Thesis, unpublished. 



14 



